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Abstract. On the occasion of the 40-th anniversary of IHES I present a few scientific 
reminiscences: most of my scientific life has been marked by my visits and I run through 
£^ ' them concluding with the analysis of a problem that originated during my last visit. The 

problem is to develop a convergent perturbative algorithm for the construction of the 
5^ ' "Eliasson's potential" for the stable and unstable manifolds of an invariant torus: and to 

Oh, study its properties. A brief review follows. 
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31. Reminiscences. 



In 1966 I came to IHES as a young "professor": I was in fact just a Ph.D. student, by 
all standards. On the recommendation of Sergio Doplicher I was taken in the group of 

■ Ph.D.'s of Daniel Kastlcr. Of course I was fully aware of "being out of place" : no time, 
however, was spent on this fact and soon I was eagerly working with Salvador Miracle- 
Sole on problems posed by David Ruelle. He went away for a few weeks and we met 
all possible difficulties so that we decided to write down the details of a proof that the 
problem given to us was not soluble: to realize, when the proof was complete, that we had 
in fact solved it. This was extending to many body interactions the Kirkwood-Salsburg 

<3\ ■ equations (also known as "cluster expansion") for lattice system: it was a tremendous 

encouragement. It was for me the opening of a new world; that of perturbation theory: 
suddenly we had in our hands a most powerful instrument. We began using it to prove, 
with Dcrcck Robinson, that the Ising model had no phase transitions at high temperature, 
uniformly in the density. 1 At about the same time Jean Lascoux arrived with several 
papers from the Russian school; we saw that our work made them transparent and 

■ we could explain them around: in particular the work of Minlos and Sinai on phase 
coexistence appeared as a great achievement and as proof of the flexibility of perturbation 
theory which could even yield the analysis of the deep two phase region. I kept thinking to 
the matter until, quite a few years later, I could understand the theory of the fluctuations 
of the interface in the 2-dimensional Ising model at low temperature. 

Leaving IHES in 1968 (I mean after may) I had met so many people (for instance Joel 
Lebowitz) and absorbed so many ideas and techniques that I was mature enough to pose 
problems without the help of Ruelle or of my collaborator and close friend Miracle-Sole. 
Saying farewell to Ruelle I saw on his blackboard (that he kept as black and powderless as 
a blackboard can possibly be until, much later, he gave it up, inexplicably) the Navier- 
Stokes equations: that he commented by stating that he was abandoning statistical 
mechanics for fluid mechanics "that seemed more interesting" . 

I brought along my copy of the fluids of Landau-Lifschitz and a few days later, in the US, 
I started studying it as I was sure I would soon need it. Even though I was discouraged by 
the difficulties (very transparent in spite of the book style which presents every problem 
as simple and as completely solved) I kept following Ruelle's work getting updates on 
the occasions of my frequent visits to IHES. Until, in 1973, I suddenly realized and was 
fascinated by his proposal of the existence of a probability distribution associated with 
chaotic motions in a conceptual generality comparable to the one, familiar in equilibrium 
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statistical mechanics, of the Boltzmann-Gibbs distribution: this was, in my view, a 
development much more significant than the (timely and necessary) critique of Landau's 
Ptolemaic theory of fluid mechanics 2 stemming out of the Ruelle-Takens paper. 
Ruelle's viewpoint on strange attractors made them simple objects via the Markov par- 
titions and Sinai's theory of Anosov systems. But I could not go further: nevertheless the 
problem ("how to obtain some directly observable consequences" from Ruelle's principle) 
remained hunting me and I thought that it was the right approach and kept lecturing 
on the subject every year since while attempting even some numerical experiments 3 but 
mostly following other people experiments. 

In the meantime I tried to understand Quantum Field Theory (I had cultivatied the 
feeling of a deep connection between renormalization theory and the Kirkwood Salsburg 
equations since the work with Miracle-Sole): in one of my visits to IHES I had met 
Francesco Guerra and in a few words, explaining his own basic work on Nelson's approach, 
he managed to make suddenly clear what I considered until then impenetrable (the works 
of Nelson and of Glimm-Jaffe and renormalization theory). So I started thinking to 
the matter and, also after a memorable lecture of Eckmann on the use of the cluster 
expansion in field theory, I was able to see the connection between the renormalization 
group of Wilson, constructive field theory and the cluster expansion developing a new 
interpretation of the mathematical works done until then in 2 and 3 dimensional field 
theory. I exposed the results at IHES with several people around (there by chance, 
because Gelfand and Sinai were visiting also) and I will never forget Niko Kuipcr calling 
me in the Director's office and telling me that my talk had been successful: the first time 
I had been in that office was when Leon Motchane, a few months after my arrival and a 
few days after completing the first paper with Miracle-Sole (before the KS-equations), 
told me that I was given a 15% raise in my salary. 

After the work in field theory I kept visiting regularly IHES where I was attracted also 
by the new young member Jurg Frolich: the policy of IHES started by Motchane of 
hiring young people at the highest rank ( "member" ) , in spite of the obvious risks, is in 
my view what made the place so exceptional and so interesting and productive. But I 
was never able to collaborate with Frolich. Instead I worked with Henri Epstein in spite 
of the age difference: together with Pierre Collet we studied smooth conjugacy between 
flows on a surface of constant negative curvature stimulated by opinions by other visitors 
that what we wanted to do was "proved" to be impossible. The effervescent atmosphere 
around Dennis Sullivan always frightened me: his synthetic and powerful approach was 
in a way opposite to my nature. But of course I was among those who were impressed 
and attracted by his view on the theory of interval maps and Feigenbaum's theory and 
tried (without success) to imitate it in the theory of the tori breakdown in conservative 
systems: a subject into which I was drawn by Joel Lebowitz' request, while we were both 
at IHES, to explain in a talk the KAM theory that I had boasted, in refereeing a paper 
for him, to have understood from Arnold's celestial mechanics paper (the seminar took 
place a few days later at the Ecole Politechnique). 

A long gap in my visits followed: which more or less coincided with a period in which I 
did not really work on new problems but looked at consequences of previous ones. I kept 
nevertheless thinking to the old questions and in my next long visit, in 1993, I was ready 
to attack a problem that E.G.D. Cohen (also a visitor at the time) proposed me to jointly 
work to find the connection between the experiments and theoretical ideas that he and 
his coworkers had been developing, in parallel with the works of Hoover and coworkers, 



Unfortunately Ptolemy has become the villain of science: I do not share at all this superficial conception 
and I consider his work as great as possible. Hence a Ptolemaic theory is very respectable in itself and 
one needs work to criticize it, if at all possible. 

Learning programming on an archcological computer located in a little room of IHES, which had a 
marvelous object, a plotter, linked to it: if Oscar Lanford, then a member of IHES, considered numerical 
investigation of chaotic motions worth of devoting time to it then I should know too. 
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on nonequilibrium statistical mechanics. Two exciting years followed: in retrospect all 
that could have been done in one afternoon; but for us it was very hard work and at 
the end of the two years we could make a precise proposal (the "chaotic hypothesis" ) for 
the application of Ruelle's principle to some simple but concrete dynamical problems. 
My last visit to IHES was in 1997: I went there with bellicose projects to continue 
working on nonequilibrium statistical mechanics. But as usual I was taken away by other 
intervening projects 4 which ended having to do with the theory of homoclinic orbits and 
the transversality of their intersections: for this reason, although this is not a problem of 
the same size of the previously mentioned ones, I will devote to it the technical part of 
this note. I am sure that, as all problems that I started at IHES, this too will be fruitful 
and linked tightly with the previous works on the cluster expansion. 
I also must remember here M.me Annie Rolland-Motchane: she, the first general secre- 
tary to IHES, clearly shared the merit of conceiving the structure of IHES and of making 
it work. 

It is a structure that during the last 33 years always made me feel that IHES was an ideal 
place as a source of inspiration and of systematic work and it has been for me a privilege 
to be able to work there. I took this literally, as privileges must be earned: I was so 
absorbed in my work, taking up most of the nights (the day being reserved to wandering 
around trying to get the most out of the people present, with a few trips to the small 
library just to remind myself that even though it was small it was still completely full of 
things that I ignored) : to the point that I could always manage avoiding the frustration 
due to my feeling of being there out of place among colleagues far more knowledgeable 
than I could possibly hope to be. 

§2. Homoclinic intersections in Hamiltonian systems. A "field theoretic" 
approach. 

Consider a Hamiltonian system: 

I 2 

H = A+ — + J g 2 (cos^ - 1) + ef(a,ip) (2.1) 

where u = {lo u uj 2 ) e R 2 , A = (A 1 ,A 2 ) £ R 2 , a = (a u a 2 ) e T 2 , I 6 R, tp 6 T 1 
and f(a,ip) is an even analytic function of the angles a,ip which is a trigonometric 
polynomial in <p; its analyticity domain is supposed to be |Imaj| < £. Therefore (2.1) 
represents a quasi-periodically forced pendulum. 

What follows can be easily extended to the system obtained from (2.1) by adding to it 
A, A) where J is a positive diagonal matrix. The case (2.1) is the limiting case 
as J — > +oo of the latter extension. The more general model is called anisochronous or 
Thirring 's model as its peculiar properties were pointed out by Thirring, [T] . 

Supposing that uj verifies a Diophantine property \u ■ v_ \ > C 1 1 i£ | 1 for all v_ G Z 2 , 
integer components nonzero vectors, it follows that for e small enough and for each A 
there is an invariant torus T(A ) which has equations parameterized by ip £ T 2 of the 
form A = A + K{±), a = ± + MV0> 1 = = S('jp). The motion on T( A ) 

is ip — > tp + tot and the average of H_ over ip (hence on time) vanishes. 

Hence we have a natural parameterization of the invariant tori by the (time) average 
value A of their action variables. In the more general Thirring's model there is also 
a relation between the rotation vector u on an invariant torus and the (time) average 
action A of the motion on it: namely u = J -1 A (this is a remarkable property which 
prompted the appellative of "twistless" for such invariant tori). 

Let E(A ) be the energy of the motions on the invariant torus T(A ). The torus 
1~(A ) is unstable and its stable and unstable manifolds W U (A ), W s ( A ) will consist 
of points with the same energy E(A ); the intersection of W a (A ), a = u,s, with the 
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plane ip = 7r and with the energy surface H = E(A ) consists of points that can denoted 
as: X(a) = (ir, a,I a (q_), A a (aj), for a = u,s and a G T 2 . In general we denote a 
point in phase space X = (<p, a, I, A). 

Therefore on the 4-dimensional "section" ip = ir,H = E(A ) the manifolds have di- 
mension 2. The vector Q(a) = A s (a) — A u (a) defines the splitting. It has been 
remarked by Eliasson that the vector Q(a) is a gradient: this means that there is a 
function <E>( a ) such that Q ( a ) = d a &(a ). I will call 3>(a) the Eliasson's potential 

Hence, since a varies on T 2 , there will be at least one point a where Q (a ) = 0. 
The corresponding point X(a ) will be, of course, homoclinic, i.e. common to both 
manifolds. In the case at hand the parity properties of / imply that a = is one such 
point, see for instance [G3]. 

It follows from the general theory of the tori T(A ), due to Graff, which is very close 
to KAM theory (particularly easy in the case (2.1), see for instance [G2]; see §5 of [CG] 
for the general anisochronous case), that the function <&(«) is analytic in e for e small, 
and such are the functions A a (a ), I a (a), for a = u,s. 

Therefore we must be able to find the coefficients of their power series expansion in e: 
in fact the expansion of A a (a ), I a (a) has been thoroughly discussed in all details in 
[G3], in the case of a polynomial /. It would be easy to derive from it the expansion for 
$(a). 

However here, §3, I will give a self contained exposition which leads to the expansion 
for &(a) borrowing from [G3] only a few algebraic identities that it would be pointless 
to prove again. And I will dedicate §4 to a brief review of the results that are known 
(to me) on the matter or that can be easily derived from existing papers (which usually 
do not deal with Eliasson's potential but rather with its gradient). A question will be 
formulated at the end of this work. 

The expansion below is the best tool, to my knowledge, to achieve a unified proof of 
various theorems dealing with the homoclinic splitting: the latter will be defined here 
as the Hessian determinant A of 3>(a) evaluated at the homoclinic point at a = 
(in our case). Evaluation of the latter determinant was begun by Melnikov who gave a 
complete solution (in spite of claims in other directions) in the "simple" case in which 
u is regarded as fixed. The quasi-periodic case was treated in the spirit of the present 
work, i.e. with attention to the Arnold's diffusion problem, in [HM] (togheter with the 
important realization of the usefulness of improper integrations analysis). 

Considerable interest has been dedicated to the problem by various authors. The results 
are not easily comparable as every author seems to give his own definition of splitting: 
the most remarkable results have been developed in the russian school approach (based 
on the key works of Melnikov, Neishtadt, Lazutkin, [Ge]). And sometimes it has even 
been difficult to realize that some papers were just plainly incorrect (like the result in 
§10 of [CG]). 

The interest of the above definition is its direct relation with the problem of showing 
the existence of heteroclinic strings, i.e. sequences of tori T(A i ) such that W u (A i ) n 
W 8 (A i+1 ) 7^ 0. Given a curve I — > A(£), with T(A(£)) being equienergetic tori 
(i.e. such that w ■ A(£) = const in the simple case (2.1), as one realizes from sym- 
metry considerations), and calling A a ( a , t), a — u, s, the equations of the manifolds of 
T(A(t)), one finds such a string if one can show that A u (a;£) — A s (a;£') — admits 
a solution for £' close enough to £. Therefore one has to apply the implicit functions 
theorem around a = and £' = £: the determinant of the Jacobian matrix controlling 
this problem is, clearly, the splitting defined above. Hence proving that the splitting is 
nonzero implies the existence of heteroclinic strings. 

In the anisochronous case not all the A {£) are necessarily the average actions of an 
invariant torus (because of the resonances always present when the rotation vectors do 
depend on the actions) so that it becomes important to measure the size of the splitting 
compared to the size of the "gaps" on the curve parameterized by £. Hence we need to 
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know quite well the dependence on u of the splitting and of the gaps and their relative 
sizes, see [GGM2]. 

I know of no paper in which the above definition of splitting is used in the case of quasi 
periodic forcing models, other than [CG] and the later [G3], [GGM1], [GGM2], [GGM3] 
or the related [BCG]: I will not discuss the papers using other definitions. 

§3. Feynman's graphs for Eliasson's potential. 

We set J = 1 for simplicity. Let a = u, s and X a (0) = (ir, a,I a (a), A a (a)); and let 
X a (t) = (<p a (a, t), cx + ut,I a (a,t), A a ( a , t)) be the solution of the equations of motion 
for (2.1): tp = I, a = a, i = -d v f a ((p) - ed,pfi(a,<p), A = -£0« /i ( a , <p) where 
/o = g 2 cos ip and /i = f(a , ip). In the case e = the stable and unstable manifolds of 
the torus with (average) action A coincide and the parametric equations of their stable 
and unstable manifolds are at ip — tt simply X = (tp = it, a, I = —2g, A = Ao) and the 
motion are X (a \t) = (<p°(t) = 4arctg e -s*, a + wt,I{t) = -gy/2(l - cos ip°(t)), A = 

Ao)- 

Hence if X a (a,t) = ( « , t) + eX a -^ ( a , t) + e 2 X a -<- 2 \ a , t) + ... we can immedi- 
ately write the equations for X^ = (tp a ^ k \ 0,I a ^A a - (k) ), k > 1, where the angle 
components vanish identically because the equations for a are trivially solved by the 
order solution. The latter is a property of the isochrony of (2.1) and it does not hold in 
the anisochronous case, which is however equally easy to treat, see [G3]. The equations 
are, dropping for simplicity the label a: 

( J (*> ) = ( g 2 cos p° J ) ' ( fw ) + ( F™ ) = L (t) ' ( J <*> ) + ( 4 fe) ) 

A {k) = £ (fc) (3.1) 

and the functions F , F_ are deduced from the equations of motion; recalling that jo = 
g 2 cosp and /i = f(a,ip): 

=i:~d^fo(p o (t)) y: n^ j) w+ 

h=2 ' k 1 + ... + k h = k j=i 

k h 

+ E-^ +1 /i(^ + ^^°w) E n^'w (3.2) 



h=l k 1 + ...+k h = k-i j = i 

kj>l 



k 



£ (k) =J2-^d^ fl (a + ut^m II 

h=l ' fci + ...+fc h =*-i j=l 

kj >i 

It is convenient to rewrite the above expressions in a more synthetic and symmetric form: 
*o W =E E ~d^f s (a + u,t,p°(t)) £ fl^Ht) 



h\ 

5=0,l/i=2-5 ki+...+k h =k-s j=l 

k s >i 

F (k) =J2 E -hd«%f s (a + u;t,<p (t)) £ f[p^\t) 

5=0,l/i=2-5 ' k 1 + ...+k h =k-s j—i 

where of course in the second relation only the 5 = 1 term does not vanish. 
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(3.3) 



Hence if W(t) = ( w ° a ^J w oi(t) j ig the solution of thc 2 x 2-matrix equation W = 



wio(t) loii(i) 

L(t)W, W(0) = 1 we find: (J') =^)(K 1 '(0) + />(t)-> ( itfV) ) ^ " 
The matrix W(t) can be easily explicitly computed: the matrix elements are holomorphic 
for .g|Imt| < -| and woo,w w tend to zero as t — > ±00 as e _fir '* while the other column 
elements tend to 00 as e 9 '*'. We shall only need the matrix elements woo, woi which have 
a simple pole at ±i^. It is w 0Q (t) = co * hgt and w i(t) = j^( c Jl gt + s'mh gt). 

Since we need only ip^ , h = l,...,k — 1, to evaluate we spell out only the 

expressions of (p( k \t) and A^{t): 



ip {k) (t)=w i(t)(lW(0)+ / w Q0 (T)F^(T)dT) - WQ0 (t) / w q1 (t)F™ (r)dr 

V t 7o 7 7o (3.4) 

A (fe) (*) = 4W(0)| / F {k \ T )d T 
Jo 

We want to determine, for each choice of a, the initial data /^(O), A^ k \0) so that thc 
motion is asymptotically quasi periodic, because we want to impose that the motion tends 
to the invariant torus T(A ). This means that the initial data must be determined in a 
different way depending on whether we impose this condition as t — > +00 or as t — > — 00: 
in the first instance we determine the stable manifold and in the second the unstable 
one. 

This condition is imposed simply by requiring A ^ (0) = — J^°° F_ ^ (r)dr where a = + 
if we impose the condition at t = +00 or a = — if we impose the condition at t = —00 and, 
likewise, /( fc )(0) = — / CT °° woai^F^^dr provided the integrals converge. Therefore 
(3.4) become: 

V {k \t) = I (woi(t)woo(T)-Woo(t)woi(T))F^ k) (T)dT- 

J (TOO 

-tuoo(t) / ^oi(r)^ (fc) (r)dT (3.5) 

J (TOO 

A {k) (t) = f F {k) { T )dT 

J (TOO 

and we have a simple recursion relation provided the integrals converge. 
The above expressions however may involve non convergent integrals and eventually 
they really do, in general. This has the consequence that (for instance) it is by no means 
true that A s '^ k \i) t ^ +oo > as the (3.5) would imply if thc integrals were proper. Of 
course there will be no ambiguity about the meaning of such improper integrals. The 
meaning is uniquely determined simply by the requirement that the asymptotic form of 
(3.4) be a quasi periodic function. 

Since the functions in the integrands can always be written as series of functions of the 
form a x lyGt ?} e -garh e iw-vt f or some x = Q,\,j,h,v with a — signr, it is easy to see that 
the rule for the evaluation of the integrals of such function is simply that of introducing 
a cut off factor e~ RaT with Re R large enough so that all the integrals converge, then one 
performs the (now convergent) integrals and then one takes the residue at R = of the 
result: see §3 of [G3] for a more detailed discussion. To derive this rule one shoould try 
a few simple cases (like evaluating the second order explitcitly). 

Improper integrals are very familiar in perturbation theory of quantum fields where 
they are normally introduced to obtain compact and systematic representations of the 
coefficients of perturbation expansions. Typically a Feynman diagram value is given by 
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an improper integral: the algorithm is so familiar that it has become usual not to even 
mention which are the rules for the evaluation of such integrals. 

Since the rules for the evaluation of the above imoroper integrals are discussed in detail 
in [G3] I shall not dwell on them and, instead, I proceed immediately to use the improper 
integrals in the same way they are used in quantum field theory: i.e. to find a simple 
diagrammatic representation of the iterative scheme described above. It is remarkable 
that such a scheme was found by Eliasson in his breakthrough theory of the KAM series, 
[E] , without any reference to field theory: he has independently developed a diagrammatic 
representation of the KAM series. 

We represent tpW (t) as: 

and, with the same "logic" : 



£(fc)( t ) 




Fig.l 



Fig.2 



where it is kj = k — S v , see (3.3); the label S v can be or 1: the first drawing 
represents the term with S = S v in the expression for F_( k > in (3.3), and the second 
drawing represents the contribution to Fg fe ' with 5 = S v . 

The node v represents —d^ +1 fs v times in the second graph and — d a dh,f\ times 
in the first. Because of the <9„ derivative we can imagine that in the first graph thelabel 
S v on the node v is constrained to be 1 . 

We can in the same way represent ip^ (t) and A ^ (t) : we can in fact change the labels t 
on the lines merging into the node v into labels t and interpret the node v as representing 
an integration operation over the time r; one get in this way the following graphs: 



A {k \t) 




(to) 



^W(t) = 




Fig.3 



>(k h ) \k h ) 

The node v with the label S v , which we noted that it must be 1 in the first drawing and 
that can be either or 1 in the second, has to be thought as representing the operations 
acting on a generic function F: 



l a F{t) - f F(t) dr (3.6) 

J (TOO 

O a F{t) = f (tUoi(t)tUoo(T) - W 0(t)w l(T))F(T)dT - W Q0 (t) f W Q1 (T)F(T)dT 
Jaoo V ' JO 

where a = + if we study the stable manifold and a = — if we study the unstable one. 
In this way the graphs of Fig.3 represent respectively the values: 



3 = 1 



(t), ^( - d h ^h[r] J] ^[r]) (t))x (3.7) 



7 



where an argument in square brackets means a dummy integration variable, inserted 
just as a reminder of the integration operation involved; furthermore fs[t] abbreviates 
f S (a + ut,ip°(t)). 

Clearly j( fe ), A (k) can be expressed simply by summing over the labels kj and S the 
values of the graphs in Fig. 3: the summations should run over the same ranges appearing 
in (3.2), i.e. h between 2 — 5 and k, and kj > 1 such that ^ . kj — k — 5 and 5 = 0, 1). If 
we study the stable manifold we must take a = + and if we study the unstable one me 
must take a = - and jW, A {k) become respectively I s + k \ A s ^ or I u ^ k \ A u ' {k) . 

It is now immediate to iterate the above representation; one simply recalls that each 
symbol: 

represents ip^ kj ^(t) and that (3.7) is multi linear in the ip( ki \t). This leads to representing 
A^ k \t) as sum of values of graphs $ of the form: 



root 




Fig.5 



A graph i? with p Vo = 2,p Vl = 2,p.u 2 = 3,p V3 = 2,p Vi = 2 and fe = 12, and some labels. The lines 
length is drawn of arbitrary size. The nodes labels 8 V are indicated only for two nodes. The lines are 
imagined oriented towards the root and each line A carries also a (not marked) label t v , if v is the node 
to which the line leads; the root line carries the label t but its "free" extreme, that we call the "root", 
is not regarded as a node. 

The meaning of the graph is recursive: all nodes v, see Fig.5, represent O a operations 
except the "first" node vq which instead represents a 1 operation; the extreme of inte- 
gration is +oo if we study the stable manifold and — oo if we stufy the unstable one. 
Furthermore each node represents a factor — ^d^ v+1 fs v [t„] if p v is the number of lines 
merging into v, except the "first" node vo which represents — da 9^ fs v [t v ] instead. The 
product Yl v i s the "combinatorial factor" for the node v 

The lines merging into a node are regarded as distinct: i.e. we imagine that they are 
labeled from 1 to p v , but we identify two graphs that can be overlapped by permuting 
suitably and independently the lines merging into the nodes. 

It is more convenient to think that all the lines are numbered from 1 to m, if the graph 
has m lines, still identifying graphs that can be overlapped under the above permutation 
operation (including the line numbers). In this way a graph with m lines will have 
a combinatorial factor simply equal to provided we define 1 instead of J[ v j-j the 
combinatorial factor of each node: we shall take the latter numbering option. Hence 
in Fig.5 one has to think that each line carries also a number label although the line 
numbers, distinguishing the lines, are not shown. 

The endnodes Vi should carry a (ki) label: but clearly unless fc, = 1 they would represent 
a (p^^ which could be further expanded; hence the graphs in Fig.5 should have the labels 
(ki) with ki — 1: this however carries no information and the labels are not drawn. The 
interpretation of the endnodes is easily seen that has to be: O t7 (—d lp fs v .)(T v ') if v\vi is 
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the line linking the cndnode i>i to the rest of the graph. An exception is the trivial case 
of the graph with only one line and one node: this represents I(—dafs vo )(t) and it will 
be called the Melnikov's graph. 

In this way we have a natural decomposition of A^ a ^ k \a,t) as a sum of values of 
graphs. It is now easy to represent the power series expansion of the trajectories on the 
manifolds W a (T( A )): one simply collects all graphs with labels S v with J2 V d v — k > 1 
(they can have at most 2k lines, if one looks at the restrictions on the labels) and adds 
up their "values" obtaining the coefficient A a,t " k \i). The O a and T a operations involve 
integrals with coo as an extreme and one has, obviously, to choose a = + if a = s and 
a = — if a = u. 

Since all the integration operations O or X are, in general, improper we see the con- 
venience of the graphical representation and its analogy with the Feynman graphs of 
quantum field theory: in fact this is more than an analogy as the above graphs can be 
regarded as the Feynman graphs of a suitable field theory: see [GGMO] for the discussion 
of a similar case (i.e. the KAM theory representation as a quantum field theory). 

An essential feature is missing: namely the graphs have no loops (they are in fact 
tree graphs). This major simplification is compensated by the major difficulty that the 
number of lines per node is unbounded (i.e. a field theory that generated the graphs 
would have to be "non polynomial"). 

Noting that the value of each graph is a function of a we now have to check that each 
A a ' {k) (0) has the form ^4 a ' (fe) (0) - 8^ {k) - 

For this purpose we consider graphs like Fig. 5 but with the root branch deleted keeping 
however a mark on the first node v$ to remember that the line has been taken away. We 
call such a graph a rootless graph. 

It is convenient to define the value Val CT (•&) of such rootless trees: its is defined as before 
but the marked node now represents the operation I a ^(F) = J° drF(r) with a = + 
for the analysis of the stable manifold and a = — for the unstable, and the function 
d p "°fSv H (keeping in mind that the marked node must have 5 Vo = 1, by construction). 

The key remark is now the identity ("Chierchia's root identity", see [G3]): 

1a,o(F O a (G)) =l„s)(GO a (F)) (3.8) 

which is an algebraic identity as our improper integrals only involve functions F, G linear 
combinations of "monomials" of the form a x IfMlie-s^'V-'- * f or SO me \ = 0, h, v_ 
with a = signr, see above, for which both sides of (3.7) can be explicitly and easily 
evaluated. 

This identity can be used to relate the values of different graphs: it means that the 
values of two rootless tress differing only because the mark is on different nodes and 
otherwise superposable are identical: this can be seen easily by successive applications 
of the identity (3.7), see [G3]. 

Therefore if we define: 

1 oo 

$ CT ' (fc) (a) = -^VaU(0), $ a (a) - e fe $"' (fe) ( a ) (3.9) 
■d fe-i 

we see that the gradient with respect to a of <& CT >( fe )(a) is precisely A a ^ k \a). And the 
splitting Q (a) is the gradient of $(a) = $ + (a) — 3>~(a). 

One can get directly a graphical representation of <&( fe ) as: 

/ — oo 
dtWal CT(t) (tf) (3.10) 
■■oo 
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where a(t) = sign (t) and Wal (T ( t ) (i?) is just the integrand in the X CTi o integral with respect 
to the first node variable r = t Vo appearing in the evaluation of Val CT ("&) . This concludes 
the construction of Eliasson's potential. 

§4. Properties of the potential. 

Many properties of the gradient Q (a) = da&((x) have been studied in [G3], [GGM1], 
[GGM2], [GGM3]: they are immediately translated into properties of the potential $, ci- 
ther by integration or by following the proofs of the corresponding statements for Q (a). 
We just summarize them: 

(1) If u is fixed then, generically, the first order dominates: 

I— oo 

*(a) = e/ drf(a + ut,<p°(t)) + 0(e 2 ) (4.1) 

this is the well known Melnikov's result. We shall say that there is "dominance of 
Mclnikov's term" for some quantity every time that the lowest order perturbative term 
gives the dominant asymptotic behavior for it in a given limiting situation. Hence in 
(4.1) domination refers to e — > 0. 

In the "one dimensional" case not explicitly treated above, but much easier, of a periodic 
forcing in which there is only one angle a and one action A, Melnikov's domination 
remains true even if the parameter g becomes small provided the e is chosen of the form 
Hg q for some q > (proportional to the degree No of / as a trigonometric polynomial in 
ip) and \/i\ small enough. 

This is somewhat nontrivial: in [G3] there is a proof based on the above formalism; 
other proofs are available as the problem is classical. The nontriviality is due to the 
necessity of showing the existence of suitable cancellations that eliminate values of graphs 
contributing to 3>(a ) higher order "corrections" (corresponding to special graphs) which 
are individually present and, in fact, larger than the first order contributions. 

(2) The next case to study is the same case of g small but with the Hamiltonian (2.1) 
(i.e. quasiperiodically, rather than periodically, and rapidly forced): let g 2 — rj and rj < 1 
be a parameter that we want to consider near 0. In this case, too, convergence requires 
that e — jirj q for some q > (proportional to the degree Nq of / as a trigonometric 
polynomial in ip) and small enough. 

The problem is discussed already in [G3] and, following it, we consider the graphs 
$ that contribute to and at each node we decompose fs into Fourier harmonics 
fs(QL,<f) !iL^~) el ~'~- This leads to considering new graphs $ in which at each 

node v a label" v v is added signifying that in the evaluation of the graph value the 
functions fs v {a,Lp) are replaced by f6 v ,v{ i p)e. 1 — '-- Of course at the end we shall have 
to sum over all the "momentum labels" v_ v G Z 2 . We call the contribution to $^ fe - ) 
from one such more decorated graph. Then from §8 of [G3] one sees that: 

v> v (4.2) 

<(b V ^) k 

where g = rji the sum over y_' runs over the nonzero values of the sums of subsets of 
v 1 , . . . , v_ k . The constants B, D are bounded by an inverse power of ij and p > is 
constant (depending on the degree of / in ip) ; the constants 6, q can be bounded only 
in terms the maximum of |/| in a strip Imaj|, |Imy>| < £ on which the maximum is 
finite. The first property follows from the analysis in §8 of [G3]; the second is symply the 
statement that the stable and unstable manifolds are analytic function of e with radius 
of convergence proportional to rj q for some q (essentially a result of Graff, see §5 of [CG]). 
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(3) A consequence of (4.2) is that can be represented as, see [GGM3] for details 
on the corresponding statement for the gradient of $: 

V_ 

M il =y,U, n (1 - cos eft (4.3) 

n J+Qc cosngt 

\d h „D„(a)\< V - p »C h \v\ h , Ph ,q h >0 

which allows us to say, very easily, that the Eliasson's potential is "in some sense" 
dominated by Melnikov's value at least in the special cases in which f„ >n are positive 
and "as large as possible" , i.e. f„ >n = ce _re l-l<5„ ; i fore > 0, and lo has good Diophantine 

def 

properties, e.g. if uj\/uj2 is the golden mean (here \v\ = \v\_\ + 

In the latter instance one verifies that, for all a, \D iL (a,e)\ < r\~ q M„ for all 1 1£ | < n^ 1 
which togheter with the analyticity in e of 3>(a) allows disregarding the contributions 
to $(a) from the z/'s exceeding r/^ 1 . The properties of the golden mean allow us 
immediately to see that in the sum only one pair ± v_ dominates at a — : it is the pair 
M-q = (fk, —fk+i) h fj is the Fibonacci's sequence and k such that n\v \ H — ^rl w • v_ \ 

is minimum; apart from exceptional intervals of values of r\ in correspondence of which 
there may be two pairs (or more) (see §2,6 of [DGJS]). The domination persists for all 
the a's such that sin t^ • a\ > b where b > is any prefixed constant (the smaller b 
the smaller has e to be to insure dominance). 

Also the gradient of $( a ), and in fact any derivative of $ is dominated by the Melnikov's 
term, by the same type of argument. But this is somewhat trivial: the real question is, 
in view of the remarks in §2 about the possible applications to heteroclinic strings and to 
Arnold's diffusion, whether the homoclinic splitting is dominated by Melnikov's integral. 
This seems to be, in the generality considered here, still an open problem. The reason is 
very simple; from (4.3) one easily deduces that: 

A = -e 2 e-(l^l+l^-<i'l)f'r i ((z, A v'fMvMv, +ed„ t „,) (4.4) 

hence one realizes that the term that should be leading, v = ±v' = ±v , is missing in 
the lowest order part. Therefore the main contribution comes, or may come, from the 
remainder d^,^ 1 on which we have little information besides the above bounds (which 
would be plenty if the Melnikov's main term did not vanish). Curiously the above 
exceptional cases, i.e. when the value of n is taken along a sequence rjj — ► such that 
for each j there are two minimizing vectors v and v_' , can be, instead, easily solved 
because (v A v[,) 2 >lasno two Fibonacci's vectors can be parallel. 
In the literature there are various claims about "proofs" of dominance of Melnikov's 
contribution to the splitting: they however seem to be always proofs of the "easy part" 
namely of the dominance of the Melnikov's term in some components of the splitting 
vector Q (a) (implied by the above analysis). 

The only known case of generic dominance of the Mclnikov term for the splitting is the 
one discussed in [GGM1], see (5) below. And its analysis is already far more subtle than 
the above. 

(4) In general the estimates (4.2), called "quasi flat" in [G3] are optimal (sec [GGM4]): 
hence one cannot hope to have bounds on the Fourier transform of $ of the form, for 
some r > 0: \Q v \ < const r]~ r e~*s\ — '-le~ K l-L Such estimates are called "exponentially 
small" and, occasionally, have been claimed to be possible. 

(5) The above results are very easy compared to the ones that can be obtained by taking 
g 1 fixed and w = {rj^uo^^^uo 1 ), discussed in [GGM1] and called the three time scales 



11 



problem, because the system has three time scales of orders respectively r]~? , . In 
this case we consider the values of rj for which uj verifies a Diophantinc property of the 
form | u ■ v_ | > if\ v_ |~* with some 7, r > and we take e equal to a suitably large power 
of 77 so that the small divisors problems can be overcome and the invariant tori do exists. 
The quasi flat estimates hold (for small rj) but they do not imply that the matrix 
da,a^(a)\a = o has three matrix elements of size exponentially small as rj — > 0. In 
fact all the four matrix elements are of the order of a power of m this is so in spite of the 
fact that the Melnikov term M(a) generates a contribution to the 2x2 splitting matrix 
with three exponentially small entries. 

In other words neither $ nor the Hessian matrix d a , a $| a — are dominated by the the 
Melnikov's "first order" contribution. Nevertheless the Melnikov 's contribution to the 
Hessian gives the leading term in the limit rj — ► 0! (generically in the perturbation). 

Of course if the above mentioned exponential estimates could be correct this would 
follow immediately from them: but they are not valid (as they would imply the wrong 
statement that the splitting matrix has 3 exponentially small entries) and the result holds 
only because remarkable cancellations take place. Hence, contrary to what is sometimes 
stated, the above case requires a delicate analysis, compared to the one in [G3] which 
solves easily the problems (l)-j-(4) above at least as far as the domination of the first order 
in the derivatives of the Eliasson's function (hence the splitting vector) is concerned. In 
particular this means that $ is not a good measurement of the splitting. 

It is in the theory of this "three time scales problem" that the analogy with field theory 
and renormalization theory turns out to be particularly useful and the methods charac- 
teristic of such theories apply very well and turn into a rather simple matter the check 
of the infinitely many identities that are necessary in order that all terms in the Hessian 
that dominate the Melnikov contribution cancel each other leaving out just the Melnikov 
contribution as the leading one as rj — > 0. 

(6) The just described graphical technique seems not only very well suited for the ques- 
tions analyzed or mentioned above but it seems quite promising also with respect to the 
solution of one of the main standing problems, namely: what is the asymptotic behavior 
of the splitting as g — > and w fixed? The case mentioned in (3) above requires that all 
the Fourier components of the perturbation do not vanish: a finer analysis shows that 
this can be somewhat weakened but not to the extent of allowing polynomial perturba- 
tions. Hence such cases seem to have a rather limited interest: they appear in fact too 
special. But even so the only thing we know is the Melnikov's dominance in Eliasson's 
potential and in its derivatives. 

On the other hand the theory discussed in [G3], [GGM3], and in §3 suggests the fol- 
lowing question. First of all let us define an extension of Melnikov's function to higher 
order. We simply consider the function $°(a) which is obtained from the diagrammatic 
representation (3.9) but replacing the operators O associated with the nodes of t by the 
operator. 

1 /"* 

Oo{F)(t) = -Y J / (w m (t)w m (T)-w m (t)w m {T))F{ T )dT (4.5) 

then, supposing w with golden rotation number (or any number with very good Dio- 
phantine nature): 

Conjecture: In model (2.1) and assuming that g — rp , will the Hessian of $>°(a ) give 
the leading asymptotics as rj — > of the splitting determinant at a = "generically" in 
f ? 

Here generic means both genericity in the space of trigonometric polynomial pcrtuba- 
tions of fixed degree (arbitrary) and in the space of the analytic perturbations, possibly 
with the constraint that the perturbation is of positive or negative type. However we 



12 



require that the perturbation be polynomial in the ip variable, see (9) below. As far as 
I know there is no proof even of the convergence of the scries defining $° (which is well 
defined only as a formal series and which may have to be regarded as an asymptotic 
series, see [G3], [GGM1]). 

The conjecture can be extended to the case of three time scales considered in (5): in 
that case it is affirmatively answered in [GGM], where, however, one also sees that the 
Eliasson's potential and its derivatives is not dominated by the first order. It is only the 
splitting determinant that is dominated by the first order: not surprisingly as this is the 
only quantity among the ones discussed which has a direct physical meaning. 

The conjecture could be strengthened by adding, for instance, that <i> can be replaced 
by the function F° obtained from <I> by developing in powers of e its Fourier coefficients 
<I>° and retaining only the lowest non vanishing order $° of each Fourier cocfficcnts to 

form the Fourier transform of <£>°. In this stronger form it becomes, in the assumptions of 
(2) above (fast forcing and "maximal size of the Fourier coefficients of the perturbation) , 
simply the statement that the splitting determinant can be computed by the first order 
Mclnikov's integral: an open problem (as mentioned above). However in this form the 
conjecture is not really stronger than above because using $° instead of $ amounts to 
saying that the edv_, u' in (4.4) has the form e{v_ A z^') 2c C v <- 

Clearly in order that the answer to the question be affirmative one has to show the 
existence of suitable cancellations: I have checked that they are indeed present at the 
order beyond the lowest (since the lowest order for the homoclinic determinant is the 
second, this means that the answer is affirmative to third order). The check requires 
using the results in [GGM1], which may already imply a positive answer to all orders. 

de f 

Denoting O the operator 0(F)(t) = O a{t) {F)(t) (with a{t) = sign (t)) one remarks that 
in all the expressions involved in the graphs evaluations one always really uses O: then 
it is useful to note the (algebraic) relation between the operator O and Oq- 



and, as it is clear from [GGM1], the G, Go factors play the role of "counterterms" in 
the field theory interpretation of the diagrammatic expansion of Hence the above 
question suggests that the leading behavior of the splitting determinant is due to graphs 
without counterterm contributions. 5 Here the "counterterms" contain non analytic func- 
tions and they are responsible for the impossibility of exponentially small estimates. A 
positive answer to the above question would state that they only give rise to subleading 
contributions to the splitting. 

An explicit expression for the value contributing to <i> can be found in [GGM1]: see 
(6.2), for the isochronous case (2.1), and see the paragraph preceding (7.4) for the 
anisochronous case. 

(7) The above theory can be immediately extended to anisochronous cases: one just 
has to consider a few new types of graphs, [G3], that contribute to the splitting vector 
Q (a) and to the splitting potential $(a ). 

(8) Most of the considerations above do not really use that the dimension of the quasi 
periodic motion is 2: if it is suppoo=scd larger it is however difficult to see what will 



Called in field theory "most divergent" graphs: rather improper an expression because in any reasonable 
field theory there should be no divergences at all; as it is the case in the theories that have been actually 
shown to exist on a mathematical basis 



OF(t) =O F{t) + \w 01 (t)\ G (F) + w 00 {t) G(F) 
°o F (t) = \ [ ( w oi(tKo(r) - w 00 (t)w 01 (T)) F{t) dr 



(4.6) 
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be the leading behavior of the splitting. One reason is that even the analysis of the 
Mclnikov term is itself a quite difficult task: Diophantine approximation theory is in a 
very rudimentary stage if the dimension of the quasi periodic motion is > 3. 

A glimpse of the difficulties that one should expect to meet is given by the three time 
scales problem (5) above. In this problem we can think that the slow frequency of order 
i]2 is in fact obtained because the perturbation by a three dimensional quasi periodic 
motion with three fast frequences wi,W2,^3 of order contains an almost resonant 
harmonic v_ such that v\ui\ + ^2^2 = 0{r]^). One would then naively think that in 
this case the homoclinic splitting can become "large" because we can have u ■ v_ small 
of order rp with not too large v\,vi. But this is illusory precisely because from the 
results of the case (5) one sees that, although we could expect a large splitting vector 
and matrix, its Hessian at the homoclinic point will be exponentially small as r\ — > 0. 
Therefore in the three dimensional case we should expect that the resonances do not 
enhance the splitting: they can make large the spliting matrix but not its determinant! 
This remark also explains why the problem (5) above is so unexpectedly difficult to 
analyze (see [GGM1]). 

(9) Finally there seems to be no reason whatsoever for having a small homoclinic splitting 
when the perturbation is not a polynomial (but "just" analytic) in the <p variable, not 
even when the rotation vector w is very fast. 
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